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Abstract 

Background: Disruptive selection has been documented in a growing number of natural populations. Yet, its 
prevalence within individual systems remains unclear. Furthermore, few studies have sought to identify the 
ecological factors that promote disruptive selection in the wild. To address these issues, we surveyed 15 
populations of Mexican spadefoot toad tadpoles, Spea multiplicata, and measured the prevalence of disruptive 
selection acting on resource-use phenotypes. We also evaluated the relationship between the strength of disruptive 
selection and the intensity of intraspecific competition — an ecological agent hypothesized to be an important 
driver of disruptive selection. 

Results: Disruptive selection was the predominant mode of quadratic selection across all populations. However, a 
directional component of selection favoring an extreme ecomorph — a distinctive carnivore morph — was also 
common. Disruptive selection was strongest in populations experiencing the most intense intraspecific competition, 
whereas stabilizing selection was only found in populations experiencing relatively weak intraspecific competition. 

Conclusions: Disruptive selection can be common in natural populations. Intraspecific competition for resources 
may be a key driver of such selection. 

Keywords: Competition, Disruptive selection, Resource polymorphism, Selection differential, Spadefoot toad, 
Quadratic selection 



Background 

Disruptive selection occurs in a population when two or 
more modal phenotypes have higher fitness than the 
intermediate phenotypes between them [1]. Disruptive 
selection has long been viewed as important in main- 
taining and increasing variation within natural popula- 
tions [1-3]; favoring the evolution of alternative 
phenotypes [4,5] and sexual dimorphism [6-8]; and even 
initiating speciation [6,8-12]. Nevertheless, compared to 
the other main modes of selection — directional selection 
and stabilizing selection — disruptive selection has trad- 
itionally received much less attention. 

Despite this relative lack of attention given to disrup- 
tive selection, recent meta-analyses suggest that disrup- 
tive selection may be at least as common as stabilizing 
selection [13,14]. Indeed, an increasing number of 
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studies have documented disruptive selection in natural 
populations [15-24]. Yet, although these data imply that 
disruptive selection may be widespread, we still do not 
know how prevalent it can be within individual systems 
where it has been found to occur [22]. 

Additionally, the actual causes of disruptive selection 
have been relatively understudied. Longstanding theory 
suggests that disruptive selection arises from negative 
frequency-dependent interactions, such as those stem- 
ming from intraspecific competition, predation, mutual- 
ism, and parasitism [10,25-28]. To illustrate how 
intraspecific competition generates disruptive selection 
[3,8,10,29,30], consider first that — owing to functional 
trade-offs [17-20,31-35] — individuals with certain pheno- 
types are generally better adapted than are individuals 
with other phenotypes at utilizing specific resource types 
[36]. Yet, as any one of these modal phenotypes becomes 
common, it will tend to suffer from resource depletion, 
as these individuals compete more against themselves 
than against other resource-use phenotypes in the same 
population [18,37]. Consequently, the less common 
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modal phenotype(s) will be favored, because they will ex- 
perience reduced competition for their resources 
[10,35,38]. In this way, intraspecific competition acts as 
an agent of frequency dependent disruptive selection, 
which favors two or more modal resource-use pheno- 
types in the same population [35,38,39]. 

Empirical studies largely support this theory. Indeed, 
intraspecific competition has been shown to cause dis- 
ruptive selection in natural populations of three-spine 
sticklebacks [22,35], spadefoot toad tadpoles [18,21], and 
Eurasian perch [23], and it has been implicated as caus- 
ing disruptive selection in various other systems [15- 
17,19,20,40]. Moreover, several studies have shown that 
competition generates negative frequency-dependence 
among different resource-use phenotypes [31,41-44], 
which (as noted above) is a hallmark of competitively 
mediated disruptive selection [38]. 

Nevertheless, relatively few comprehensive surveys of 
natural populations have sought to identify the eco- 
logical conditions that are associated with disruptive se- 
lection [22,23]. For example, the relationship between 
spatial variation in the strength of intraspecific resource 
competition and the occurrence and intensity of disrup- 
tive selection has been little examined in un- 
manipulated wild populations. This relationship has 
likely been difficult to evaluate because doing so requires 
evaluating selection in replicated populations, some of 
which experience weak resource competition and some 
of which experience strong resource competition. 

We sought to fill these gaps in our knowledge concerning 
the prevalence and correlates of disruptive selection in nat- 
ural populations. We specifically tested two predictions: 
first, that disruptive selection may be common within cer- 
tain systems, and second, that disruptive selection would be 
more intense under conditions of greater resource competi- 
tion. We tested these two predictions in natural populations 
of Mexican spadefoot toads {Spea multiplicata). 

Spadefoot toad tadpoles are well-suited for such stud- 
ies, because they express a remarkable range of trophic 
phenotypes in the wild [45,46]. The extremes of this 
variation are represented by two ecomorphs that com- 
prise a resource polymorphism: an "omnivore" morph — 
a round-bodied tadpole with a long intestine, small jaw 
muscles, numerous rows of keratinized denticles, and 
smooth keratinized mouthparts that feeds primarily on 
the pond bottom, and a "carnivore" morph — a narrow- 
bodied tadpole with a short intestine, greatly enlarged 
jaw muscles, few rows of keratinized denticles, and 
notched, serrated keratinized mouthparts that feeds 
mostly in the water column [18,47,48]. Omnivores are 
generalists, which feed mostly on microscopic detritus, 
algae, and small crustaceans, whereas carnivores are 
specialists, which feed mostly on anostracan fairy 
shrimp [37,49]. 



Although carnivore development is induced by con- 
sumption of fairy shrimp [50], heritable variation for 
morph development exists within natural populations 
[48,51]. The degree to which a population expresses this 
resource polymorphism varies across ponds, and is, in 
part, associated with variation in conspecific density and 
ecological opportunity (i.e., the presence of underutilized, 
accessible resources). Specifically, bimodality in trophic 
phenotypes is greatest in ponds where conspecific density 
and ecological opportunity are highest [46] . 

Previous research suggests that intermediate phenotypes 
are disfavored by disruptive selection in this system 
[18,21,46]. In particular, compared to tadpoles with inter- 
mediate phenotypes, omnivores and carnivores are larger, 
more developmentally advanced, and more likely to survive 
to metamorphosis [18,21]. Furthermore, previous experi- 
ments have shown that this disruptive selection reflects 
negative frequency-dependent interactions driven by eco- 
logical specialization and resource competition [43,46]. 

However, it is unclear how prevalent disruptive selection 
is among populations of S. multiplicata or how ecological 
variation impacts the mode and strength of selection in 
this system. We therefore addressed these issues by meas- 
uring the mode and magnitude of phenotypic selection, as 
well as the strength of intraspecific competition within 
natural populations. We focused on populations of Mexi- 
can spadefoot toad tadpoles (S. multiplicata) in the San 
Simon Valley of southeastern Arizona and southwestern 
New Mexico, USA (Figure 1). These populations are ideal 
for addressing the above issues because the 5. multiplicata 
in these populations express a resource polymorphism 
that is thought to have arisen from disruptive selection 
[18,21,46]. Moreover, previous research suggests that dis- 
ruptive selection may be widespread in these populations 
[18,21]. Finally, the 5. multiplicata in these populations 
experience a wide array of ecological conditions over a 
small geographic area [21,52,53]. 

We made a priori predictions based on theory as well 
as the prior work on this system (see above). Specifically, 
we predicted that disruptive selection on tadpole trophic 
morphology would be widespread in our surveyed popu- 
lations. We also predicted that this disruptive selection 
would be strongest in ponds where intraspecific compe- 
tition is the most intense (as measured by conspecific 
density and per capita resource abundance). 

Results 

Evaluating the prevalence of disruptive selection 

Disruptive selection was the predominant mode of quad- 
ratic selection on tadpole trophic morphology (Figure 1). 
Indeed, disruptive selection, identified by a significandy 
positive y (Table 1) and a fitness minimum (Figure 2), oc- 
curred at least once in 73% of the unique ponds sampled 
(11 of 15 unique ponds) and in 59% of our total 
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Figure 1 Localities. Map of study area illustrating locations of ponds sampled and the form of quadratic selection in each pond. Symbols for 
ponds sampled in multiple years are divided to show the form of quadratic selection in each year. The numbers beside each symbol section 
correspond to each collection's pond ID referenced as "Map ID" in Tables. 



collections (Table 1; 13 of 22 collections). In contrast, sta- 
bilizing selection, identified by a significantly negative y 
(Table 1) and a fitness maximum (Figure 2), occurred in 
13% of the unique ponds sampled (2 of 15 unique ponds) 
and in 14% of our total pond collections (Figure 1, Table 1; 
3 of 22 collections). Additionally, a significant directional 
component of selection favoring carnivore-like phenotypes 
was present in 80% of the unique ponds sampled (Table 1; 
12 of 15 unique ponds) and 82% of total pond collections 
(Table 1; 18 of 22 collections). 

Evaluating whether competition predicts the strength of 
disruptive selection 

As predicted, disruptive selection was more intense in 
ponds with greater intraspecific competition (i.e., lower 



per capita resource density; Figure 3A, F 1>2 o = 9.861, 
P = .005, median regression coefficient from bootstrap- 
ping =-.016, 95% CI: -.023 - -.011, higher conspecific 
density; Figure 3B, F h20 = 10.088, ^=.004, median re- 
gression coefficient from bootstrapping = .023 95% CI: 
.016 - .030). 

Discussion 

Interest in disruptive selection arises from its potential 
role in maintaining and accentuating variation within 
populations (including polymorphism) and in promoting 
speciation [2]. Despite this longstanding interest, the 
prevalence of disruptive selection in natural populations, 
and therefore the relative importance of disruptive selec- 
tion in generating biological diversity, remains unclear. 
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Table 1 The mode and strength of selection on trophic morphology in natural ponds 



Pond 


Map ID 


Year 


B/v 


Selection differentials 


SE/2SE 


f 


P 


AZ0602 


1 


2006 


p 


.022 


.006 


3.678 


.0003 








Y 


.013 


.012 


1.181 


.240 


AZ0603 


2 


2006 


P 


.020 


.005 


3.667 


.0003 








Y 


.019 


.011 


1.639 


.103 


AZ0604 


3 


2006 


P 


.014 


.008 


1.622 


.111 








Y 


.039 


.013 


2.874 


.006 


AZ0605 


4 


2006 


P 


.031 


.009 


3.487 


.0006 








Y 


.126 


.014 


8.258 


<.0001 


AZ0606 


5 


2006 


P 


.003 


.002 


1.423 


.159 








Y 


-.001 


.004 


-.170 


.866 


AZ0607 


6 


2006 


P 


-.017 


.007 


-2.354 


.020 








Y 


.038 


.012 


2.920 


.004 


NM0608 


/ 


2006 


P 


.009 


.004 


2.253 


.025 








Y 


.024 


.006 


3.596 


.0004 


AZ0710 


8 


2007 


P 


.026 


.006 


3.857 


.0002 








Y 


.043 


.008 


4.552 


<.0001 


AZ0711 


9 


2007 


P 


.076 


.009 


8.262 


<.0001 








Y 


-.060 


.018 


-3.113 


.002 


AZ0706 


10 


2007 


P 


.006 


.003 


1.67 


.097 








Y 


.005 


.005 


.936 


.351 


AZ0801 


1 1 


2008 


P 


.030 


.006 


4.417 


<.0001 








Y 


-.036 


.008 


-4.618 


<.0001 


AZ0810 


12 


2008 


P 


.042 


.006 


6.841 


<.0001 








Y 


.075 


.012 


6.467 


<.0001 


AZ081 1 


13 


2008 


P 


.008 


.003 


2.318 


.021 








Y 


.024 


.004 


5.741 


<.0001 


AZ0816 


14 


2008 


P 


.029 


.004 


6.781 


<.0001 








Y 


.068 


.006 


11.412 


<.0001 


AZ0809 


15 


2008 


P 


.017 


.006 


2.796 


<.007 








Y 


.009 


.008 


1.224 


.226 


AZ0802 


16 


2008 


P 


.005 


.001 


2.714 


.007 








Y 


.011 


.002 


4.377 


<.0001 


AZ0812 


17 


2008 


P 


.009 


.003 


2.699 


.007 








Y 


.029 


.006 


5.557 


<.0001 


AZ0813 


18 


2008 


P 


.017 


.005 


3.095 


.002 








Y 


.028 


.01 


3.110 


.002 


NM0810 


19 


2008 


P 


.006 


.004 


1.518 


.131 








Y 


.033 


.008 


3.812 


.0001 


AZ0903 


20 


2009 


P 


.020 


.012 


1.638 


.105 








Y 


.001 


.020 


.089 


.929 
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Table 1 The mode and strength of selection on trophic morphology in natural ponds (Continued) 



AZ0902 


21 


2009 


P 


.017 


.002 


6.017 


<.0001 








Y 


.016 


.006 


2.692 


.007 


AZ0904 


22 


2009 


P 


.012 


.002 


5.326 


<.0001 








Y 


-.010 


.002 


-4.917 


<.0001 



The pond name, map ID corresponding to Figure 1, and year of collection are given for each population, along with the regression terms (P/y), estimated 
selection differential for each term, its standard error (SE), f-statistic and probability of rejecting the null hypothesis that the estimated differential is zero. For 
quadratic regressions, positive selection differentials signify disruptive selection and the quadratic regression coefficient is doubled to calculate the quadratic 
selection differential (y) and the associated standard error (SE) is also doubled. Bolding signifies statistical significance. 



We sought to determine the prevalence of disruptive se- 
lection among natural populations of spadefoot toads. 
We also sought to establish whether the strength of dis- 
ruptive selection was positively associated with the in- 
tensity of intraspecific competition within individual 
populations. 

Disruptive selection favoring extreme resource-use 
phenotypes, as measured here by using larval body size 
to estimate tadpole fitness (see Methods), was wide- 
spread in the spadefoot toad populations surveyed in 



this study. Indeed, we documented disruptive selection 
on tadpole trophic morphology in 11 of 15 populations 
surveyed (Figure 1). These results, together with those of 
previous studies [18,21,46], strongly support the role of 
disruptive selection in the evolution of resource poly- 
morphism in Spea tadpoles [46], and add to the evidence 
that disruptive selection can be widespread within cer- 
tain systems [22]. 

Although disruptive selection was prevalent in our sur- 
veyed populations, we also found considerable variation 




Morphological index 

Figure 2 Fitness functions. Cubic-splines of relative fitness (measured by tadpole body size) on a composite shape variable of trophic 
morphology. The cubic spline (solid line) is bracketed by 95% confidence intervals (dashed lines). Individual panel legends correspond to the 
populations' map ID followed by the pond ID. An asterisk indicates a significant fit of a quadratic regression, as well as a fitness minimum or 
maximum. 
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Figure 3 Intensity of disruptive selection as function of the intensity of competition. The relationship between quadratic selection and two 
measures of intraspecific competition. The intensity of disruptive selection on trophic morphology was greater in ponds with (A) lower per capita 
resource density and (B) higher conspecific density. Standardized quadratic selection differentials above the dotted line are positive and indicate 
disruptive selection, while those below the line are negative and indicate stabilizing selection. Significant selection differentials (black), non- 
significant selection differentials (grey) and significant selection differentials with no fitness minimum or maximum in the range of the data 
(open) are all shown. Solid and dashed lines show respectively, the median and 95% confidence intervals of the regression coefficient obtained 
from bootstrapping. Each point in the analysis was weighted by the populations' sample size. 



in both the mode and magnitude of quadratic selection 
(Figures 1, 2; Table 1). This variation correlated predict- 
ably with the intensity of intraspecific competition 
within individual populations. Specifically, disruptive se- 
lection was strongest in ponds with the highest conspe- 
cific density and lowest per capita resource density 
(Figure 3 A, B). In contrast, stabilizing selection favoring 
intermediate resource-use phenotypes was found in 
ponds where intraspecific competition was weak 
(Figure 3A, B). These results therefore confirm theory 
suggesting that intraspecific competition can drive dis- 
ruptive selection (see Introduction). Furthermore, our 
results suggest that disruptive selection is likely to occur 
in habitats relatively free of predators, heterospecific 
competitors, or other ecological factors that depress 
population size and that therefore potentially weaken 
intraspecific competition. 

Interestingly, the two extreme ecomorphs do not ap- 
pear to be equally favored by disruptive selection. Previ- 
ous studies found that carnivores achieve greater 
survival and are generally (although not always) larger at 
metamorphosis than either intermediates or omnivores 
[18,21]. We similarly found that a directional component 
to selection favoring carnivores may be widespread 
(Figure 2; Table 1). Yet, if one extreme morph (the car- 
nivore morph) has higher fitness than the other ex- 
treme morph (the omnivore morph), how is resource 
polymorphism maintained in this system? 

Several factors examined in previous studies may help 
maintain both omnivores and carnivores within the 



same population. First, as carnivores become increas- 
ingly common in any given pond, and as competition 
among carnivores for the limited shrimp resource 
thereby becomes more intense, negative frequency- 
dependent selection favors omnivores [43]. In support of 
these earlier findings, we found evidence in this study 
for pure disruptive selection (with no directional compo- 
nent) in two ponds, and disruptive selection with a dir- 
ectional component favoring omnivores in another 
(Figure 2; Table 1). Moreover, although carnivores bene- 
fit by gaining access to the more profitable shrimp re- 
source (which could explain why they often achieve 
larger body size), additional experiments have shown 
that carnivores also experience greater competition with 
other carnivores than omnivores do with other omni- 
vores [37]. In this way, negative frequency-dependent se- 
lection maintains both morphs within the same 
population [43]. Second, fitness trade-offs are associated 
with each of these phenotypic alternatives. Compared to 
carnivores, omnivores invest more into abdominal fat 
bodies, which increases post-metamorphic resistance to 
starvation [43]. Thus, although carnivores may often 
have higher survival before metamorphosis (the life stage 
at which we estimated selection in the present study), 
omnivores may have higher survival immediately after 
metamorphosis, possibly balancing their lower pre- 
metamorphic survival [43]. Generally, frequency- 
dependent selection, coupled with fitness trade-offs, 
likely contribute to the evolutionary maintenance of 
many resource polymorphisms [4]. 
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A deeper problem requiring explanation, however, is 
why intermediate phenotypes persist in the face of wide- 
spread disruptive selection against them (Figure 1). Al- 
though morph determination is environmentally 
triggered — carnivore development is induced when a 
tadpole consumes fairy shrimp [50]— considerable herit- 
able variation in the response to this cue exists in nat- 
ural populations [48,51], suggesting that the propensity 
to produce intermediate phenotypes could, in principle, 
be eliminated by natural selection. A possible resolution 
to this problem is that sexual selection may lead to ran- 
dom mating in regard to larval phenotype, recreating 
tadpoles likely to develop as intermediate phenotypes in 
each generation. Specifically, in our populations, female 
S. multiplicata prefer to mate with males in good condi- 
tion, irrespective of the females' own condition [54,55]. 
Therefore, female S. multiplicata will tend to mate with 
males of good condition, potentially regardless of either 
the male's larval phenotype or the female's own larval 
phenotype (e.g., a good-condition adult male might have 
been either an omnivore or a carnivore as a tadpole). 
This directional selection on male condition (exerted by 
female preferences for mates) might thereby oppose dis- 
ruptive selection acting on the larvae. Indeed, directional 
sexual selection may often oppose disruptive natural se- 
lection in many populations [56], which may explain 
why an increasing number of cases have been reported 
in which assortative mating by resource-use phenotype 
was expected but not found [57,58]. 

Spatial and temporal variation in the mode of selection 
may also help explain the persistence of intermediate 
larval phenotypes. For example, stabilizing selection 
favoring individuals with intermediate phenotypes was 
detected in some ponds (Figure 1, Figure 2), specifically 
those lacking strong intraspecific competition 
(Figure 3A, 3B). Stabilizing selection was also detected 
in ponds in which S. multiplicata co-occurs with a con- 
gener, Spea bombifrons [21]. Such mixed-species ponds 
may be within a few kilometers of pure S. multiplicata 
ponds in the San Simon Valley of southeastern Arizona 
[53]. Thus, S. bombifrons migrants may occasionally 
colonize nearby pure S. multiplicata ponds, thereby 
changing the selective regime within the pond to one in 
which intermediate phenotypes are favored. Indeed, S. 
bombifrons occurred along with S. multiplicata within 
some of the ponds sampled for this study as recently as 
30 years ago [59]. Thus, because the mode and direction 
of selection is spatially and temporally variable, inter- 
mediate phenotypes may persist in the population. 

Finally, limits in the ability of tadpoles to respond ap- 
propriately through phenotypic plasticity to variation in 
both fairy shrimp density [50] and their competitive en- 
vironment (i.e., conspecific density and morph fre- 
quency) may explain the persistence of intermediate 



phenotypes in natural ponds. For instance, spatial and 
temporal variation in the density of fairy shrimp within a 
pond could give some individuals a head start in carni- 
vore expression early in development. Once tadpoles in- 
crease their mobility as they grow and develop, tadpoles 
encountering fairy shrimp later may begin to develop as 
carnivores, but these late-developing carnivores might 
end up being outcompeted by earlier-developing carni- 
vores. If these individuals are deprived of fairy shrimp, 
their carnivore-like features begin to regress, and they 
assume an intermediate phenotype [43]. Generally, inter- 
mediate phenotypes may be difficult to eradicate via dis- 
ruptive selection when limits exist in the ability of 
individuals to assess and respond adaptively to their 
competitive environment though phenotypic plasticity. 
Further studies are needed to test these ideas. 

Conclusions 

Although disruptive selection has long been viewed as a 
potential generator of biological diversity, it has received 
relatively little attention, and, consequently, little is 
known regarding its prevalence or causes within specific 
systems. Using spadefoot toad tadpoles as a model sys- 
tem, we found that disruptive selection can be prevalent 
in the wild, and that its occurrence is predictably asso- 
ciated with a specific ecological factor: strong intraspeci- 
fic competition. Given that intraspecific competition is 
common in many natural populations — and frequently 
strong [60,61] — our results therefore imply that disrup- 
tive selection may be a more important force contribut- 
ing to the origin and maintenance of biological diversity 
than is currently appreciated. 

Methods 

All procedures were carried out in compliance with the 
Institutional Animal Care and Use Committee at the 
University of North Carolina under application 06- 
047.0-A. Field collections were conducted under New 
Mexico collecting permit 1857 and Arizona collecting 
permit SP604848. 

Field surveys 

We collected S. multiplicata tadpoles during summers 
2006-2009 from 15 ephemeral ponds in the San Simon 
Valley (Figure 1). In all 15 ponds, S. multiplicata was the 
only species of Spea present. Six ponds were sampled 
across multiple years for a total of 22 collections 
(Figure 1; Table 2). We sampled each pond on a single 
day 16-20 days after breeding had occurred (Spea breed 
only a single time within a pond each season, shortly 
after the pond fills with rainwater). Within each pond, 
we sampled tadpoles from five randomly selected sites 
throughout the pond using a hand-held dip net (median 
sample size = 130 tadpoles). We sacrificed the tadpoles 
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Table 2 Summary of tadpole collections and ecological parameters in natural ponds 



Pond 


Map ID 


Year 


N 


Tadpole density 


Fairy shrimp density 


Cover (%) 


AZ0602 


1 


2006 


93 


2 


3 


.8 


AZ0603 


2 


2006 


124 


1 


2 


.75 


AZ0604 


3 


2006 


50 


2 


3 


.5 


AZ0605 


4 


2006 


176 


3 


2 


.6 


AZ0606 


5 


2006 


94 


2 


1 


.9 


AZ0607 


6 


2006 


102 


2 


1 


.4 


NM0608 


/ 


2006 


165 


3 


2 


1 


AZ0710 


8 


2007 


78 


3 


1 


.9 


AZ071 1 


9 


2007 


99 


2 


3 


1 


AZ0706 


10 


2007 


125 


1 


2 


./ 


AZ0801 


11 


2008 


99 


1 


2 


1 


AZ0810 


12 


2008 


213 


3 


2 


.8 


AZ081 1 


13 


2008 


181 


3 


3 


.5 


AZ0816 


14 


2008 


297 


3 


2 


1 


AZUoUy 


1 C 

1 J 


zUUo 


jy 




3 


i 


AZ0802 


16 


2008 


150 


2 


2 


./ 


AZ0812 


17 


2008 


188 


2 


2 


.6 


AZ0813 


18 


2008 


135 


2 


2 


1 


NM0810 


19 


2008 


169 


2 


2 


.8 


AZ0903 


20 


2009 


78 


3 


I 


.75 


AZ0902 


21 


2009 


211 


3 


3 


1 


AZ0904 


22 


2009 


192 


1 


1 


1 



For each population the pond name, map ID corresponding to Figure 1, sampling year, tadpole sample size (N), tadpole density, fairy shrimp density, and 
percentage of vegetative cover (Cover) are shown. We assigned numerical values to our estimates of tadpole and shrimp abundance such that "high" = 3, 
"medium" = 2 and "low" = 1 . 



immediately after collection by immersion in a 0.1% 
aqueous solution of tricane methanesulfonate (MS 222) 
and preserved them in 95% ethanol. We used this ran- 
dom sampling technique to estimate the density of S. 
multiplicata tadpoles in each pond as 'high| 'moderate' 
and 'low' [46]. Additionally, we determined the range of 
available resources in each pond by estimating abun- 
dances of fairy shrimp and detritus, which are the two 
main resources on which Spea tadpoles feed. We esti- 
mated fairy shrimp abundance by sweeping a net 
throughout each pond and categorizing shrimp densities 
as "high," "moderate," and, "low". These subjective esti- 
mates were corroborated by previously published inten- 
sive, quantitative sampling [50,53]. We assessed the 
availability of detritus by estimating the percent vegeta- 
tive cover in a twenty meter radius around each pond's 
circumference and categorized each pond as having ei- 
ther "high" (67%-100% cover), "moderate" (34%-66% 
cover) or, "low" (0%-33% cover) detritus (ponds with 
more vegetation tend to have more detritus; see [53]). 

We calculated an estimate of per capita resource abun- 
dance for each pond by dividing the sum of fairy shrimp 
and detritus abundance ("high" = 3, "medium" = 2 and 



"low" = 1) by S. multiplicata tadpole density ("high" = 3, 
"medium" = 2 and "low" = 1). We subsequently used this 
score to test our prediction that ponds with high levels 
of intraspecific competition experienced the most in- 
tense disruptive selection. 

Evaluating the prevalence of disruptive selection 

To evaluate the prevalence of disruptive selection on 
trophic morphology across ponds, we first calculated a 
composite index of trophic morphology, separately for 
each pond, following previously described methods 
[18,21,46]. Briefly, we began by measuring each tadpole's 
snout-vent-length (SVL) using hand-held digital calipers. 
For each tadpole, we additionally measured the width of 
the orbitohyoideus (OH) muscle and characterized the 
shape of each tadpole's keratinized mouthparts (MP) on 
an ordinal scale from one (most omnivore-like) to five 
(most carnivore-like), and counted the number of den- 
ticle rows (DR) surrounding the keratinized mouthparts. 
We standardized OH for body size (SVL) by regressing 
In (i.e., natural log) OH on In SVL and used the resulting 
residuals for the subsequent analyses (these residuals 
were distributed normally). We then combined the MP 
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Table 3 Principal component analysis of trophic morphology 



Pond Map ID PC MP eigenvector loading residual OH eigenvector loading DR eigenvector loading % variance explained 

AZ0602 1 1 .675 .692 -.252 49.5 

2 -.252 -.104 -.961 32.4 

3 -.692 .713 .104 18.0 
AZ0603 2 1 .707 .707 — 82.0 

2 .707 -.707 — 17.9 

AZ0604 3 1 .594 .647 -.477 62.5 

2 -.511 -.154 -.845 25.7 

3 .620 -.746 -.238 11.6 
AZ0605 4 1 .576 .644 -.502 59.9 

2 -.581 -.108 -.806 25.9 

3 .574 -.756 -.312 14.2 
AZ0606 5 1 .707 .707 — 76.5 

2 .707 -.707 — 23.5 

AZ0607 6 1 .691 -.466 -.372 52.2 

2 -.057 -.057 -.882 31.7 

3 -.720 -.882 -.287 15.7 
NM0608 7 1 .701 .690 -.176 52.3 

2 .058 -.190 -.979 32.8 

3 .710 .697 -.093 14.7 
AZ0710 8 1 .600 .591 -.538 64.0 

2 -.318 -.441 -.838 20.8 

3 -.733 .675 -.076 15.14 
AZ0711 9 1 .590 .475 -.652 64.3 

2 .530 -.837 -.130 25.7 

3 -.608 -.269 -.746 9.8 
AZ0706 10 1 .560 .608 -.561 71.0 

2 .712 -.007 .701 17.8 

3 -.423 .793 .438 11.1 
AZ0801 11 1 .534 .589 -.605 44.8 

2 .836 -.470 .279 28.7 

3 -.120 -.656 -.745 26.3 
AZ0810 12 1 .588 .547 -.594 68.0 

2 .442 -.833 -.330 18.4 

3 .676 .068 .732 13.4 
AZ0811 13 1 .597 .590 -.542 53.5 

2 -.327 -.437 -.837 24.9 

3 .731 -.678 .068 21.5 
AZ0816 14 1 .597 .551 -.581 66.3 

2 .242 -.815 -.525 19.0 

3 -.764 .172 -.621 14.5 
AZ0809 15 1 .593 .558 -.578 49.4 

2 -.211 .802 .557 26.1 
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Table 3 Principal component analysis of trophic morphology (Continued) 



3 


-.776 


.208 


-.595 


24.4 


AZ0802 16 1 


.628 


.684 


-.369 


55.9 


2 


-.429 


-.090 


-.898 


30.9 


3 


-.648 


.723 


.237 


13.1 


AZ0812 17 1 


.641 


.501 


-.580 


59.9 


2 


-.125 


.815 


.565 


25.6 


3 


.756 


-.290 


.585 


14.4 


AZ0813 18 1 


.550 


.777 


-.304 


38.7 


2 


-.607 


.122 


-.784 


36.6 


3 


-.572 


.617 


.539 


24.6 


NM0810 19 1 


.577 


.595 


-.558 


76.0 


2 


-.571 


-.192 


-.797 


14.3 


3 


.582 


-.779 


-.229 


9.6 


AZ0903 20 1 


.404 


.710 


-.576 


47.5 


2 


.814 


.006 


.579 


33.5 




-.415 


.703 


.576 


1 O.O 


AZ0902 21 1 


.591 


.556 


-.583 


76.5 


2 


.307 


-.824 


-.475 


14.0 


3 


.745 


-.101 


.658 


9.4 


AZ0904 22 1 


.557 


.588 


-.585 


49.2 


2 


.828 


-.350 


.436 


26.1 


3 


-.051 


.728 


.683 


24.5 


For each population the pond name, 


map ID corresponding to Figure 1, principal component axis (PC), trait loadings and % variance 


explained are shown. 



scores, DR counts, and residuals of In OH regressed on 
In SVL, into a single multivariate shape variable (the 
"morphological index"; see [18]) with a principal compo- 
nent analysis using the correlation matrix, and standar- 
dized to have a standard deviation of one. Our 
morphological index for each pond was the first princi- 
pal component. Larger values of PCI correspond to 
more carnivore-like tadpoles, with larger OH muscles 
(corrected for body size), fewer denticle rows, and more 
serrated, notched mouthparts. In contrast smaller values 
of PCI correspond to more omnivore-like tadpoles, with 
smaller OH muscles (corrected for body size), more den- 
ticle rows and smooth mouthparts (Table 3). 

We were unable to obtain denticle row counts for 
tadpoles from two collections (AZ0603 and AZ0606, 
see Table 3). Therefore, for these two collections we 
calculated an alternative morphological index using 
only MP scores and residual values of In OH corrected 
for In SVL [21]. We calculated a single multivariate 
shape variable by combining the MP scores and resi- 
duals OH values by again using a principal component 
analysis as described above. As before, the morpho- 
logical index for each of the two ponds was the first 
principal component. Larger values of PCI correspond 
to more carnivore-like tadpoles, with larger OH 



muscles (corrected for body size), and more serrated, 
notched mouthparts. In contrast smaller values of PCI 
correspond to more omnivore-like tadpoles, with smal- 
ler OH muscles (corrected for body size), and smooth 
mouthparts (Table 3). 

We estimated the mode and magnitude of selection 
acting on tadpole trophic morphology in each pond 
using body size as a fitness proxy [In SVL; see 21, 62]. 
Body size is positively correlated with fitness in many 
species [62], including in Spea. For example, relative to 
tadpoles within the same population that are larger, tad- 
poles that are smaller in body size have a lower probabil- 
ity of survival, both before [63] and after metamorphosis 
[64]. Moreover, smaller tadpoles tend to be less develop- 
mentally advanced, and there is a premium on rapid de- 
velopment in the ephemeral ponds in which Spea 
multiplicata typically breed [43]. Additionally, in Spea, 
adult size is positively correlated with mating success in 
males [54] and fecundity in females [63]. Finally, a previ- 
ous field experiment, in which tadpoles of different mor- 
photypes were marked and recaptured within a natural 
pond, established that the relationship between morph- 
ology and body size mirrored that between body size 
and survival [18]. In short, body size is a reliable proxy 
for fitness in this system. 
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Table 4 Mitchell-Olds and Shaw constrained regression 
tests for fitness minimum/maximum indicating quadratic 
selection 



Pond 


Map ID 




p . 

' mm 


^max 


p 

' max 


fitness 
min/max 


AZ0604 


3 


5.124 


.028 


10.522 


.002 


Y 


AZ0605 


4 


48.257 


<,0001 


80.282 


<.0001 


Y 


AZ0607 


6 


12.177 


.0007 


5.311 


.024 


Y 


NM0608 


/ 


7.810 


.005 


16.573 


<.0001 


Y 


AZ0710 


8 


1 2.492 


.0006 


30.103 


<.0001 


Y 


AZ071 1 


9 


28.468 


C0001 


1.559 


.214 


N 


AZ0801 


11 


31.686 


C0001 


8.232 


.005 


Y 


AZ0810 


12 


23.971 


<,0001 


61.126 


<.0001 


Y 


AZ081 1 


13 


21.384 


<,0001 


38.224 


<.0001 


Y 


AZ0816 


14 


85.242 


<,0001 


172.109 


<.0001 


Y 


AZ0802 


16 


24.155 


C0001 


10.773 


.001 


Y 


AZ0812 


17 


20.838 


C0001 


36.537 


<.0001 


Y 


AZ0813 


18 


5.154 


.024 


14.562 


.0002 


Y 


NM0810 


19 


10.912 


.001 


16.500 


<.0001 


Y 


AZ0903 


20 


.083 


.773 


.317 


.574 


N 


AZ0902 


21 


1.691 


.194 


14.560 


.0001 


N 


AZ0904 


22 


39.010 


C0001 


9.252 


.002 


Y 



For each population the pond name, map ID corresponding to Figure 1, test 
statistics and results evaluating if the null hypotheses that fitness minimum/ 
maximum lie at extreme minimum and maximum phenotypic values (rather 
than within the observed range of the data) can be rejected. 



To estimate the magnitude and mode of selection, we 
ran separate linear and quadratic regression (including 
both linear and quadratic terms for the latter) of relative 
fitness (In SVL/mean In SVL) onto the morphological 
index for each pond [65,66]. We obtained standardized 
selection differentials for linear and quadratic selection 
from these regressions of relative fitness on trophic 
morphology. We doubled the quadratic regression coef- 
ficients to obtain quadratic selection differentials (y) 
[66,67]. 

A significant, standardized linear selection differential 
(C) indicates that directional selection is acting on 
trophic morphology [66]. Selection on trophic morph- 
ology might be disruptive when the quadratic selection 
differential (y) is significantly positive, and stabilizing 
when it is negative [66]. However, a significant y is ne- 
cessary, but not sufficient to indicate the presence of 
quadratic selection where intermediate phenotypes are 
at the fitness maximum/minimum [68]. Therefore, for 
each pond, we fit cubic splines between the morpho- 
logical index and fitness along with 95% confidence 
intervals to visualize the selective surface [69]. Cubic 
spline analysis is less sensitive to outliers and allows esti- 
mation of a fitness function without an a priori 



assumption about the function's shape [69]. We visually 
inspected these plots to determine if there was an inter- 
mediate fitness minimum or maximum within the range 
of the data (Figure 2). For the ponds with significant 
quadratic regression differentials we also applied a con- 
strained regression method to statistically evaluate if the 
null hypotheses that fitness minimum/maximum lie at 
extreme phenotypic values (rather than within the 
observed range of the data) can be rejected (Table 4) 
[68,70]. Where the graphical and statistical methods 
conflicted (in two cases) we favored our visual evalu- 
ation of the cubic splines because the constrained re- 
gression method is sensitive to deviations from the 
assumptions of parametric statistical tests that do not 
affect the nonparametric cubic spline approach [69] . 

Evaluating whether competition predicts the strength of 
disruptive selection 

Finally, we sought to test our prediction regarding the 
ecological correlates of disruptive selection. Specifically, 
we predicted that disruptive selection would be strongest 
in ponds with intense intraspecific competition. To test 
this prediction, we fit separate linear regressions of the 
quadratic selection differentials onto two estimates of 
intraspecific competition: (1) conspecific density, and (2) 
per capita resource density. Our underlying assumption 
was that intraspecific competition would be more in- 
tense the higher the density of conspecifics and the 
lower the per capita resource density within a pond. 

We weighted the selection differentials in each regres- 
sion described above by the square root of our sample 
size for each pond. We did so because confidence in the 
estimation of both the sign and magnitude of selection 
differentials is lower for those ponds with small sample 
sizes then for those with larger sample sizes. 

To further evaluate our models we present confidence 
intervals obtained from bootstrapping for each analysis to 
account for potential pseudo-replication introduced by 
using temporal replicates from the same pond. For the re- 
gression models used in testing our two predictions, we 
sampled a single selection differential from each of the fif- 
teen unique ponds, fit the regression, and extracted the re- 
gression coefficients. We resampled the selection 
differentials and refit each model 1000 times and then 
estimated 95% confidence intervals for each regression co- 
efficient. We assessed significance by asking if the range of 
the confidence interval excluded zero. All statistical ana- 
lyses were performed using R (version 2.15.0) [71]. 
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